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Abstract 



Properties of low-variability periods in the time series are analysed. The theoretical 
approach is used to show the relationship between the multi-scaling of low- variability 
C^ ■ periods and multi-affinity of the time series. It is shown that this technically simple 

P^ \ method is capable of reveling more details about time series than the traditional 

F^ I multi-affine analysis. We have applied this scaling analysis to financial time series: 

Ph ' a number of daily currency and stock index time series. The results show a good 

^ , scaling behaviour for different model parameters. The analysis of high-frequency 

USD-EUR exchange rate data confirmed the theoretical expectations. 
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1 Introduction 



a 

O ■ Perhaps the only thing which is more complicated and unpredictable than 

a human mind, is the collective human mind. This collective human mind 
is the driving force of financial market fluctuations. The intrinsic complexity 
r> \ of the market dynamics forms the basis for the exponential growth of the 

c^ I econophysics, since the publication of seminal papers of Stanley et al [1,2,3]. 

As in every new scientiflc discipline, the Econophysics is developing rapidly 
in various directions. Research in those various flelds has contributed several 
stylized statistical properties of asset returns (which, in principle, may cease 
to be valid at a certain moment of time, because the collective human mind is 
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capable of learning), c.f. [4,5,6]. Most advancements in econophysics are owing 
to the concepts of scale-invariance, intermittency and lieteroscedasticity. So, 
viable approximations for the market movements have been either based on 
the continuous time random walk (CTRW) model [7,8,9,10,11,12], or derived 
from it, c.f. the model of fractional Brownian motion (fBm) in multifractal 
trading time (fBm-mtt) [13,14,15]. There is strong empirical evidence that 
fluctuations in financial markets possess multifractal statistics [16,17,18,19,20]. 
Multifractality has been found to be inherent both to the currency markets 
[16,18,19], and to the stock exchange markets [20]. The time series analysis is 
usually based on the closing prices of two consecutive trading days, but higher 
frequency data have been also used [21]. 

The focus of this paper is placed on the multi-scaling properties of currency 
and stock exchange time series. So, we refrain ourselves from disussing the 
advancements in understanding the other aspects of financial time series. 

It has been recently shown that in the case of intermittent time series, the 
scaling behaviour of low-variability periods can provide additional informa- 
tion, as compared to a multifractal analysis [22]. Preliminary study of the 
time series of stock indices and currency rates indicated that similarly to the 
heart rate variability data, the financial time series are also characterized by 
a multi-scaling behaviour of low- variability (i.e. "silent" or "calm") periods 
[23]. Such a behaviour has been independently verified by Kaizoji et al [24]. 

In this paper, analytic approach is used to derive the relationship between the 
multi-scaling of low-variability periods, and multi-affinity of the time series. 
Besides, a detailed analysis of the statistics of low variability periods is per- 
formed for various daily stock indices and daily currency rates. The study of 
the high-frequency data of the USDEUR exchange rate indicates that at the 
time-scale of one day, there is a cross-over between two scaling laws for the 
low-variability periods. 



2 Theoretical approach to low-variability periods 



Most modern studies of the variability of intermittent time series are based on 
(or derived from / related to) the model of multi-affine fractional Brownian 
motion (mafBm), cf. [4,13,25,26,27]. This is not surprising, because multifrac- 
tal behaviour is believed to be the most universal case of scale-invariance, c.f. 
[28] . Note that while CTRW and fBm-mtt models are somewhat more detailed 
than a data description by multifractal formalism (by introducing the concept 
of waiting- or trading time), the time series generated by both these models 
can be also characterized by the multi-affine spectra and lead to the same 
scaling of low-variablity periods as the model of mafBm. 



2.1 Multi-affine fractional Brownian motion 



One of the disadvantages of the multi-affine analysis is neglecting the long- 
term correlations in the dynamics of high-frequency variability. This claim is 
motivated as follows. Consider a multi-affine time series x{t), with a lower 
cut-off scale tq. Then, each point of the time series to is characterized by its 
own Hurst exponent h (referred to as the Lipschitz- Holder exponent); this 
exponent describes the local scaling of fluctuations, 



|x(t)-a;(to)| ~ |t-to|^ |t-to|>T-o. 



Further, the distribution of points of certain values of h is self-similar, and is 
described by a fractal dimension f{h). So, according to the mafBm model, we 
study the scaling of increments \x{t) — x{to)\ at |t — to| ^ tq, without specific 
attention to the values of them at |t — to| ~ tq. In order to shed light into 
this phenomenon, the method of scaling analysis of low- variability periods has 
been devised [22]. 

On the other hand, a low-variability period of length /j is defined as such a 
continuous time interval Ti = [ti,ti + k] {i = 1,2, . . .) that (a) 

\x{t) - {x{t))^ \<5 ioT teTi (2) 



where 5 is a threshold parameter and angular braces denote sliding average 
over a window of width r > tq (in principle, the window width can be arbitrar- 
ily large; however, the highest time resolution of the metod and widest scaling 
range is achieved when it is as small as possible, i.e. just few cut-off scales tq); 
(b) each period has maximal possible length, implying that decreasing tj or 
increasing /j would lead to violation of Eq. 2. We speak about multiscaling be- 
haviour, if the cumulative distribution function of the low-variability periods 
(the number of periods with /j > n) scales as 

R{n) = Ron~''^^'^\ (3) 



where a{6, r) is a scaling exponent and Rq is a constant. 

It should be noted that the scaling exponent a{6, r) can be considered as the 
extension of the concept of multi-affine spectrum f{h) to the highest possible 
time resolution. This is because of two circumstances. First, the effective time 
resolution of the multi-affine description by Eq. 1 is limited by the required 
scaling range |t — to| 3> tq; meanwhile, the effective time resolution of the 
description based on low-variabihty periods is the cut-off scale Tq. Second, in 



the case of infinite time resolution, tliere is one-to-one correspondence between 
tlie spectra f{h) and a{S,T), see below. 

The relationship between f{h) and a{6,T) spectra can be derived as follows. 
If a time moment t is characterized by a Lipschitz-Holder exponent h, the 
mismatch between the local average and the instant value is of the same order 
of magnitude as the variability amplitude of x{t) over the time window of 
width r: 

kW - {x{t))^ I ^ t\ (4) 



Comparing Eqns 2 and 4, we conclude that for a period Ti, all the points 
are characterized hj h < log^ 6. Similarly, the edges ti — Tq and ti + k + 
Tq of the period Tj are characterized hj h > log^ 6. Therefore, the the low- 
variability periods are the contiguous pieces which remain from the entire 
time-axis after the removal of the high-variability points with h > log^ 6. 
These high- variability points form a sparse Cantor-dust-like set, assuming that 
log^ S > ho, where Hq corresponds to the global maximum of the spectrum 
f{h) [i.e. /(/io)=l]. Indeed, f{h) is the fractal dimension of the set of points 
described by the exponent h. So, the fractal dimension of the set of points 
described by exponents h > log^ 6 is found as 



. 1, if log^5 < hn 

d = sup f{h) = <^ ' Br - y ^^^ 

'»>i°sr'5 /(log^(5), if \og^5>ho 



If an interval is divided into pieces by a Cantor dust of dimension d, the 
cumulative length-distribution exponent of the pieces is also d (because the 
number of pieces larger than n is estimated as the number of boxes of size 
n required to cover the Cantor dust, according to the box-counting method). 
So, we have correspondence 

a{6,T) = f {log, 6) if \og,S>ho. (6) 



Note that in the case of simple fBm (which is studied in more details in Ref. 
[29]), there is no scaling of low- variability periods with lengths / > r. Indeed, 
for simple fBm, f{h) = 1, ii h = H [otherwise, f{h) is not defined]. So, the 
high-variability points either are absent (i.e. there is a single low-variability 
period), or they are populated quasi-homogeneously with fractal dimension 
d = 1 (then, all the low- variability periods are shorter than or of the order of 
the window width t). 



2.2 Multi-affine trading time and discussions 



The approach used for mafBin can be also apphed to the fBm-mtt model 
with x{t) = F{T{t)), where F is a fractional Brownian function of its argu- 
ment (described by the Hurst exponent H) and T{t) is a multi-affine function 
of time. Indeed, in that case Eq. 4 is substituted by 

kW - {<t))^ I ^ (A.T)^ ^ r^^. (7) 



Here, A^-T ^ r^ is the increment of the trading time T{t) and h is the local 
Lipschitz-Holder exponent for the function T(t). Accordingly, Eq. 6 is substi- 
tuted by 

a{6,r) = f{H~Hog,6) if log, 6 > Hho- (8) 



So, both in the case of the mafBm and fBm-mtt models, the scaling described 
by multifractal spectrum f{h) with h > ho can be also described by the scaling 
of low-variability periods. The approach to describing the Lipschitz-Holder 
exponent range h < ho is straightforward: instead of low-variability periods, 
high variability periods have to be used. Then, Eqns. 6 and 8 will remain 
valid, except that the inequalities will be opposite; details of this aspect will 
be analysed elsewhere. Now, we are left with two questions. 

First, the multifractal dimension f{h) cannot exceed the topological dimension 
one. Is it possible to have the scaling exponent of low-variability periods larger 
than one? The answer is yes (examples are provided below). More specifically, 
equality a > 1 assumes that the most part of the aggregated time Tagg of 
the low-variability periods corresponds to shortest periods (because Tagg = 
Jn ■ dR{n) = Ro J n'^'^dn; here, the integral vanishes at large lengths n). 
This assumes that the fractal dimension of the high variability periods is one, 
i.e. log^S < ho- In that case, the low-variability moments form a Cantor-dust 
like set of points. For real time series, the time resolution is always finite; 
then, in the regions of high "dust density" , the "dust particles" may overlap 
and form relatively long (much longer than the time resolution, but much 
shorter than the total time interval) low-variability periods. Such a clustering 
of small-/i-points is not described by the multifractal spectrum f{h). 

Second, the multifractal spectrum f{h) is a one-parameter-curve; meanwhile, 
a{6, t) is a two-parameter-curve, which is effectively reduced to one degree of 
freedom by Eq. 6 [or 8] . Is it possible to have a such a class of scale-invariant 
time series, which is more generic than the class of multi-affine functions, so 
that the exponent a can be arbitrary function of two independent variables 6 
and r. This question is left open for further studies. 



To conclude, the scaling analysis of low variability periods provides a simple 
and superior (in the sense of time resolution and applicability to the param- 
eter range with a > 1) alternative to the traditional multi-afiine analysis of 
the time series. In particular, Eq. 5 provides a way to check the validity of the 
assumption of multi-affinity, which requires a{6,T) = a{log^6). If there is no 
need for such verification, the route is even simpler, e.g. the analysis can be 
performed at a fixed window length r. All this does not mean that studying 
the scaling of low-variability periods is always better than the multi-fractal 
formalism. Instead, the two methods should be considered as complementing 
each other, because the accuracy and applicability of both methods are sensi- 
tive with respect to the width of the scaling range and specifics of the scaling 
behaviour. 



3 Properties of low-variability periods in financial time series 



Now we apply the above developed theory to financial time series. The specifics 
of long-term stock-market time series is that prices can vary by orders of 
magnitude. Therefore, it makes sense to work with a logarithmic scale, and 
substitute the definition of low variability (Eq. 2) as follows: 

|l-P(t)/(P(t))J<5 for ten (9) 



Here, unlike in the case of Eq. 2, it is assumed that average is taken not over a 
time window centered around the current moment of time t, but instead, the 
large-t-edge of the window coincides with t (because for practical applications 
in real time, the future is not known): 

{p{t))^ = -j:p{t-k). (10) 



So, a low-variability period is the period between two consecutive price move- 
ments that are exceeding the given threshold S. The length of a low variability 
period is measured in the same units as time t — in the units of the sampling 
period of the market data. The local average in Eqn 9 gives us the second 
input parameter: r that could be interpreted as a degree of "locality". With 
these two input variables we maintain the easy interpretation of the results: 
threshold value 6 is simply the percentage change of asset at time t compared 
against the average price of that asset in preceding period with length of r 
time units. Further we measure the lengths of low-variability periods and we 
count them. Then we define a cumulative distribution function R{n) that rep- 
resents the amounts of periods where low-variability lasted at least n time 
units. Finally, we fit the data against the scaling law given by Eq. 3. 



Following analysis will give answers to two essential questions: 

(1) Do the low variability periods obey power law and if so, what is the 
scaling exponent? 

(2) Assuming that there is a power law, how does the scaling exponent de- 
pend on the parameters 6 and r? 



3. 1 Determination of Power law 



We used the following data in our analysis: 



Table 1 

Data used in the analysis 

Abbrv Description 

EUR EUR / USD exchange rate 

AUD AUD / USD exchange rate 

CAD CAD / USD exchange rate 

DEM* DEM / USD exchange rate 

FRF* FRF / USD exchange rate 

GBP GBP / USD exchange rate 

JPY JPY / USD exchange rate 

DAX German equity index - DAX 

UKX UK equity index - FTSE 1 00 

NKY Japanese equity index - Nikkei 225 

FAFF MSCI Europe Austrasia Far East index 

World MSCI World index 

SPX Standard & Poors ' 500 

DJIA Dow Jones Industrial average 



Calendar period 

01/01/2001-31/12/2003 

13/12/1983 - 20/04/2004 

04/01/1971-20/04/2004 

04/01/1971-20/04/2004 

04/01/1971-20/04/2004 

04/01/1971-20/04/2004 

04/01/1971-20/04/2004 

01/10/1959-20/04/2004 

03/01/1984-20/04/2004 

05/01/1970-20/04/2004 

31/12/1969-20/04/2004 

31/12/1969-20/04/2004 

30/12/1927-20/04/2004 

03/01/1900-20/04/2004 



Frequenc 


:y # of data 


1. min 


1106031 


Daily 


5261 


Daily 


8491 


Daily 


8490 


Daily 


8491 


Daily 


8476 


Daily 


8485 


Daily 


11169 


Daily 


5134 


Daily 


8492 


Daily 


8552 


Daily 


8553 


Daily 


16767 


Daily 


26206 



*Remark: from January, 1st 1999 German Mark and French Franc are fixed 
against Euro. 

In our analysis, the scaling exponent was found by using least-squared fit for 
the whole data series. For visualizations, the dependence R{n) is plotted in log- 
log scale. Then, the scaling exponent a is equal to slope of the fit-line. Financial 
time series are single realizations of intermittently fiuctuating nonstationary 
time series, which makes calculation of exact error estimation of the scaling 
exponent impossible. However, rough estimates of the uncertainties have been 
obtained as follows. The least-squares fitted trend-line was found as described 
above, except that the slope a was not optimized, i.e. it was considered as 
a fixed parameter. Further, the sum of squared residuals r(a) was calculated 
as a function of a. The error estimate was found as e = {a' — a), where a 
is the least-squares fitted value of the slope, and a' satisfies the condition 
r{a') = 2r{a). The technique used for obtaining the scaling exponent a is 
illustrated in Fig 1. 



Above-mentioned procedure was carried out for all of the time series and 



Fig. 1. Determination of the scaling exponent a for SPX time series under the 
fohowing conditions: (5=1.6% and r=10 days 
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with different thresholds and window widths. The quality of the data fit was 
measured by the correlation coefficient R^. Total 1004 measurements were 
carried out with all the daily time series. In Fig 2, the validity of power law 
hypothesis is demonstrated. Histogram shows that in most cases, the power 
law provided a good data fit. Calculated values of a and error estimations for 
currencies and equity indices for r = 10 days are presented in tables 2 and 3 
respectively. 

It is important to outline two effects: 

• There is reverse relationship between the correlation coefficient and the 
threshold parameter 5. The reason is that the amount of large movements 
is small and not representative for statistical analysis. 

• There is also reverse relationship between the correlation coefficient and 
the window width r, due to similar reasons (longer averaging smoothes 
time series and large movements become more frequent). 

Due to these effects there is no sense to measure the scaling exponents beyond 
certain parameter range, e.g. for r = 2 days and 5 > 2%. 



Fig. 2. Histogram of R-squared coefficients based on the regression analysis of de- 
termining scaling exponent a 
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Table 2 



Scaling exponents 
Threshold S(%) 


0.10 


2.33 


0.25 


1.91 


0.40 


1.64 


0.55 


1.55 


0.70 


1.46 


0.85 


1.35 


1.00 


1.31 


1.15 


1.16 


1.30 


1.09 


1.45 


1.06 


1.60 


0.75 


1.75 


0.73 


1.90 


0.57 


2.05 


0.51 


2.20 


0.51 


2.35 


0.51 


2.50 


0.54 


2.65 


0.45 


2.80 


0.40 


2.95 


0.40 



for currency time 
CAD AUD 

+ 0.25 3.41+0.47 

+ 0.09 2.60 + 0.25 

+ 0.11 2.37 + 0.19 

+ 0.09 2.27 + 0.24 

+ 0.08 1.95 + 0.16 

+ 0.09 1.98 + 0.19 

+ 0.07 1.80 + 0.18 

+ 0.06 1.67 + 0.13 

+ 0.06 1.48 + 0.08 

+ 0.08 1.59 + 0.09 

+ 0.07 1.52 + 0.12 

+ 0.07 1.39 + 0.12 

+ 0.07 1.24 + 0.13 

+ 0.06 1.27 + 0.12 

+ 0.06 1.14 ±0.13 

+ 0.05 1.13 + 0.06 

+ 0.05 1.05 + 0.07 

+ 0.04 0.95 + 0.05 

+ 0.04 0.92 + 0.04 

+ 0.04 0.79 + 0.04 



series with r 
GBP 

1.38 + 0.07 

1.35 + 0.04 

1.51+0.05 

1.61+0.06 

1.72 + 0.08 

1.80 + 0.10 

1.71+0.08 

1.77 + 0.09 

1.77 + 0.11 

1.62 ±0.11 

1.52 ±0.10 

1.44 ±0.09 

1.36±0.11 

1.02 ±0.07 

1.14±0.09 1.19 

1.19 ±0.09 1.32 

1.01 ±0.06 

1.11 ±0.05 

0.99 ± 0.05 



10 days 
FRF 



1.38 
1.28 
1.23 
1.35 
1.35 
1.43 
1.46 
1.43 
1.38 
1.43 
1.49 
1.49 
1.51 
1.45 



1.34 
1.17 
0.97 



0.80 ±0.04 0.94 



±0.07 
±0.08 
±0.05 
±0.07 
±0.06 
±0.06 
±0.06 
±0.07 
±0.06 
±0.06 
±0.06 
±0.08 
±0.10 
±0.12 
±0.12 
±0.11 
±0.10 
±0.08 
±0.05 
±0.05 



1.88 
1.70 
1.64 
1.64 
1.69 
1.57 
1.66 
1.72 
1.65 
1.58 
1.53 
1.54 
1.43 
1.46 
1.44 
1.21 
1.10 
1.09 
0.98 
0.95 



DEM 

±0.12 
±0.11 
±0.05 
±0.07 
±0.05 
±0.07 
±0.08 
±0.08 
±0.09 
±0.07 
±0.07 
±0.08 
±0.06 
±0.08 
±0.09 
±0.07 
±0.05 
±0.06 
±0.04 
±0.05 



JPY 

1.12 ±0.06 
1.43 ±0.05 
1.36 ±0.06 
1.48 ±0.06 
1.54 ±0.05 

1.58 ±0.06 
1.60 ±0.07 
1.65 ±0.09 

1.59 ±0.08 

1.47 ±0.09 

1.48 ±0.10 
1.47 ±0.08 
1.41 ±0.08 
1.27 ±0.06 
1.24 ±0.07 
1.26 ±0.07 
1.26 ±0.07 

1.24 ±0.07 

1.25 ±0.08 
1.22 ±0.08 



3.2 Dependence of the scaling exponent on the parameters 



Two phenomena can be observed: 

(1) Exponent a and threshold 6 tend to be negatively related: larger 6 values 
correspond to lower a values. 

(2) For equities, the dependence of the exponent a on the threshold 6 is 



Table 3 

Scaling 

Threshold 

0.10 

0.25 

0.40 

0.55 

0.70 

0.85 

1.00 

1.15 

1.30 

1.45 

1.60 

1.75 

1.90 

2.05 

2.20 

2.35 

2.50 

2.65 

2.80 

2.95 



exponents for 



S{%) 



UKX 

3.87 ±0.32 
3.02 + 0.19 
3.09 + 0.30 
2.93 ± 0.34 
2.54 + 0.38 

2.44 + 0.43 
2.38 + 0.26 
2.04 + 0.20 
1.92 + 0.19 
1.76 + 0.20 

1.88 + 0.23 
1.62 + 0.12 

1.58 + 0.13 
1.70 + 0.11 

1.59 + 0.12 

1.45 + 0.11 
1.51+0.10 
1.40 + 0.08 
1.38 + 0.07 
1.35 + 0.08 



equity index 
DAX 

4.37 + 0.38 4.07 
3.87 + 0.43 3.34 
3.07 + 0.31 
3.13 + 0.35 
2.85 + 0.29 
2.72 + 0.24 
2.54 + 0.24 
2.36 + 0.23 
2.26 + 0.22 
2.21+0.25 
2.26 + 0.22 
1.97 + 0.17 
1.82 + 0.11 
1.82 + 0.12 
1.82 + 0.10 
1.78 + 0.09 
1.80 + 0.11 
1.42 + 0.08 
1.53 + 0.10 
1.53 + 0.10 



time series with 



2.64 
2.75 
2.35 
2.21 
2.20 
1.91 
1.97 
1.92 
1.82 
1.60 
1.59 
1.63 
1.38 
1.38 
1.32 
1.41 
1.34 
1.28 



NKY 
+ 0.49 
+ 0.33 
+ 0.24 
+ 0.27 
+ 0.20 
+ 0.17 
+ 0.22 
+ 0.12 
+ 0.13 
+ 0.11 
+ 0.09 
+ 0.07 
+ 0.07 
+ 0.07 
+ 0.07 
+ 0.09 
+ 0.08 
+ 0.07 
+ 0.08 
+ 0.07 



EAFE 

4.27 + 0.56 

3.16 + 0.49 
2.58 + 0.31 
2.43 + 0.22 
2.43 + 0.19 
2.36 + 0.21 
1.97 + 0.17 

2.17 + 0.17 
1.86 + 0.15 
1.96 + 0.17 
1.83 + 0.13 
1.55 + 0.11 
1.39 + 0.06 
1.35 + 0.05 
1.39 + 0.06 
1.39 + 0.06 
1.38 + 0.07 
1.31+0.08 
1.30 + 0.07 
1.24 + 0.07 



r=10 days 

World 
50 + 0.36 
98 + 0.32 
74 + 0.37 
52 + 0.26 
63 + 0.28 
34 + 0.20 
27 + 0.17 
12 + 0.15 
11+0.16 
90 + 0.12 
67 + 0.11 
54 + 0.07 
45 + 0.08 
47 + 0.09 

36 + 0.09 
18 + 0.11 

37 + 0.10 
30 + 0.06 
30 + 0.06 
09 + 0.05 



4.17 
3.32 
3.04 
2.95 
2.60 
2.39 
2.33 
2.08 
2.03 
1.97 
1.95 
1.79 
1.73 
1.65 
1.46 
1.68 
1.63 
1.50 
1.32 
1.42 



SPX 
+ 0.51 
+ 0.23 
+ 0.23 
+ 0.31 
+ 0.19 
+ 0.16 
+ 0.13 
+ 0.11 
+ 0.11 
+ 0.09 
+ 0.12 
+ 0.08 
+ 0.08 
+ 0.09 
+ 0.08 
+ 0.11 
+ 0.14 
+ 0.10 
+ 0.10 
+ 0.11 



DJIA 

4.03 + 0.57 
3.75 + 0.48 
3.09 + 0.35 
2.84 + 0.32 
2.72 ± 0.40 
2.74 + 0.26 
2.51+0.21 

2.37 + 0.15 
2.41+0.17 
2.28 + 0.15 
2.26 + 0.19 
2.19 + 0.13 
1.98 + 0.10 
1.90-0.13 
1.80 + 0.09 
1.83 + 0.11 
1.77 + 0.13 
1.77 + 0.15 

1.38 + 0.08 
1.55 + 0.13 



relatively strong, but relatively weak for currencies. 

In order to illustrate this observation, the values of a are plotted in Fig 3 
against the threshold 6 for DEM and SPX time series. 



Fig. 3. a values for DEM and SPX under thresholds 6 = 0.10% ...2.95% and r=10 
days. 
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Next we study the dependence of a on the window width r using the high 
frequency data of EUR-USD exchange rate. As mentioned above, low amount 
of statistics becomes an issue, when the threshold parameter is too high for 
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the given window width. Therefore only those values of a are plotted in Fig 4, 
where I^ > 0.7. 



Fig. 4. Scaling exponent a for EUR time series with different window widths r under 

thresholds 5= 0.10%.. .2.95%. 
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3. 3 Discussions of the results 



According to the theoretical expectations, a should be a decreasing function 
of (5, owing to Eq. 6 and to the circumstance that for h > Hq (i.e. for a < 1), 
f{h) is a decreasing function of h. This trend is expected to held for the range 
q; > 1, as well, due to simple reasoning. Increasing the threshold parameter 
6 leads to some fluctuations ceasing to be qualified as "large" , so that longer 
low-variability periods will emerge. Increasing the number of long periods, in 
its turn, corresponds to decreasing the scaling exponent a, in accordance with 
the results depicted in Fig 3. 

In order to study, how well the theoretically derived relationship Eq. 6 (or 
Eq. 8) holds for the financial time series, the scaling exponent a{6, r) has been 
plotted versus log^S in Fig. 5 (using the same data as in Fig. 4). Note that 
the data-points with r > 2880 min = 2 days describe the scaling behaviour 
at the time scale above one day; the data-points in the range r < 1 day are 
influenced significantly by the intra-day dynamics of the exchange rate. 
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Fig. 5. Scaling exponent a for EUR time series with different values of r and 5 are 
plotted versus log^<5, providing a data collapse for r > Iday, a < 1. 
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For the range of validity of Eq. 6 (a < 1), the data-points with r > 2 lay 
reasonably well into a single curve, confirming the validity of Eq. 6. The data 
of Fig. 5 indicate also that at the time-scale around one day, there is a cross- 
over between different scaling regimes: the data with r < 1 day incorporate 
a mixed dynamics (which cannot be described by mafBm model) and lay no 
longer on the same curve. At the range a > 1, there is no evident reason for 
the data to lay on a curve: all the corresponding data points are scattered in 
a form of a disperse cloud. 



4 Conclusion 



In this paper a new method for time series analysis is derived and analytically 
motivated. It is shown that for multi-affine time series, the scaling properties 
of low-variability periods are described by scaling exponent a as a function of 
threshold parameter 5 and averaging window width r. The scaling analysis of 
low variability periods offers a simple alternative to the multi-affine analysis of 
the time series, providing wider applicability range and somewhat higher time 
resolution. An open question is about the existence of a class of scale-invariant 
time series, more generic than multi-affine ones, violating Eq. 6. 

Particular emphasize is paid to the financial time series. The analysis showed a 
good scaling behaviour of the currency exchange rate and stock index data and 
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confirmed the theoretical expectations, particularly the relationship between 
the scaling exponent a and the model parameters 5 and r (c.f. Eq. 6). The 
analysis of high-frequency data of EUR-USD exchange rate indicated that 
there is a cross-over between different scaling regimes at the time scale of 
one day. It was also observed that scaling exponent values for equity time 
series tend to be more sensitive with respect to the threshold parameter than 
currency time series. 
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